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ABSTRACT 



O 

'. We study a field-theoretical analogue of the Aharonov-Bohm 

effect in two-, three- and four- dimensional Abelian Higgs models; the 
corresponding topological interaction is proportional to the linking 
number of the Abrikosov vortex and the particle world trajectories. 
We show that the Aharonov-Bohm effect gives rise to a nontrivial 
interaction of charged test particles. The numerical calculations in 
the three-dimensional model confirm this fact. 

X 

g ! 1 Introduction 

It is well known that the Abelian Higgs model in three and four dimensions has classical 
solutions called the Abrikosov-Nielsen-Olesen strings (vortices) These strings carry 
quantized magnetic flux (vorticity), and the wave function of the charged particle which 
is scattered on the string acquires an additional phase. The shift in the phase is a 
physical effect which is the field-theoretical analogue of the quantum-mechanical 
Aharonov-Bohm effect ||: strings play the role of solenoids that scatter the charged 
particles. The topological long-range Aharonov-Bohm interaction between strings and 
charged particles was discussed in the continuum theory |2| and on the lattice ||]. 

In Section 2 we show that in the lattice Abelian Higgs model with non-compact 
gauge field the Aharonov-Bohm effect gives rise to long-range Coulomb-like interaction 
between test particles. In the three-dimensional case the induced potential is confining, 
since the Coulomb interaction grows logarithmically. In Section 3 we give the results 
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of numerical calculations of the induced potential. These results confirm the existence 
of the Aharonov-Bohm effect in the Abelian Higgs model. 



2 Potential Induced by the Aharonov— Bohm Effect 

The partition function of the D-dimensional non-compact lattice Abelian Higgs mo- 
del isQ: 

+ 00 +7T 

Z= Jva Jv v £ exp{-S{A,<p,l)} , (1) 



where 



S(A, <p, I) = (3\\dA\\ 2 + 7 ||d<^ + 2nl -NA\\ 2 , (2) 

A is the non-compact gauge field, <p is the phase of the Higgs field and I is the integer- 
valued one-form. For simplicity, we consider the limit of infinite Higgs boson mass, the 
radial part |$| of the Higgs field, l^e^ being frozen. For the interaction of the Higgs 
field with the gauge field we use the Villain form of the action. 

One can rewrite Q the partition function (|XJ) as a sum over closed vortex trajectories 
(D = 3) or as a sum over closed string world sheets (D — 4), using the analogue of 
Berezinski-Kosterlitz-Thauless (BKT) transformation [0: 

ZocZ BKT = £ expf-^^'^A + m 2 )- 1 ^')}, (3) 

5*j=0 

where m 2 = N 2/ y//3 is the classical mass of the vector boson A. The closed objects 
*j defined on the dual lattice have dimension 1 and correspond to a closed vortex for 
D = 3, whereas for D = 4 they have dimension 2 and correspond to a closed string 
world sheet. These objects *j, which are topological defects on the lattice, interact 
with each other through the Yukawa forces (A + m 2 ) 

Applying the same transformation to the quantum average of the Wilson loop 
for the test particle of charge M, Wm(C) = exp{— iM(A, jc)}, we get the following 
formula 01: 



m j 



< W M (C) > n = £ exp{-47r 2 7 (*j, (A + 

*j(*c 2 )e% 

8*j=0 

~Uc (A + m 2 Y^ c ) - 2m^(* 3c , (A + m 2 )~Vj) + 2ni^E (j,j c )} . (4) 



1 We use the notation of the calculus of differential forms on the lattice B; see Appendix A. 
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The first three terms in this expression are short-range Yukawa forces between defects 
(strings or vortices) and test particles of charge M. The last long-range term has a 
topological origin: K/(j,jc) is the linking number between the world trajectories of the 
defects *j and the world trajectory of the test particle jc: 

LU,jc) = Cjc,AT 1 d*j). (5) 

In three (four) dimensions, the trajectory of the vortex (string) *j is a closed loop 
(surface) and the linking number iL is equal to the number of points at which the loop 
jc intersects the two (three) dimensional volume bounded by the loop (surface) *j. 
Equation ([5]) is the lattice analogue of the Gauss formula for the linking number. This 
topological interaction corresponds to the Aharonov-Bohm effect in field theory |2|, |j, 
the defects (vortices or strings) which carry the magnetic flux 2ir/N scatter the test 
particle of charge M. 

In the limit N 2/ y ^> /3, the partition function ([3|) becomes 



Z, 



BKT 



E 

*i(*c 2 )e^ 

<5*i=0 



exp 



Air 2 (3 



(6) 



In the corresponding 3D (4D) continuum theory the term ||*j|| 2 is proportional to the 
length (area) of the world trajectory *j; therefore, the vortices (strings) are free. In 
the said limit, the expectation value of the Wilson loop Wm(C) 



is 



< W M {C) > 



1 



E ex p{- 



N 2 BKT 

*j(*C 2 )G^ 
5*j=0 



47T 2 /3 lla , ll2 M , 



(7) 



This formula describes the Aharonov-Bohm interaction of free defects carrying mag- 
netic flux 27r/N with the test particle of charge M. 

For D = 2, the expression ([7]) can be computed exactly. In this case, the defects are 
pointlike, *j is attached to sites (*j = *j(*C2)) and the condition S*j = is satisfied 
for any *j; therefore, all the variables *j in equations (^|) and (0) are independent. The 
linking number can be written as K/(j,jc) = (j, mc)i where the two-form mc represents 
the surface spanned on the contour C. After these remarks the evaluation of Wm{C) 
becomes trivial and the result is as follows: 



< W M {C) > 



( E exp {- 



4tt 2 /3 -2 
' N 2 J 



5(C) 



V 



Eex P {-^} 

j€2Z J 



/ 



= const. exp{-a{M,N;f3) -5(C)}, (8) 
where S(C) is the area of the surface bounded by the contour C. The string tension is 
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<t(M,N;P) = -In 



re( 


'47r 2 /3 
N 2 ' 


M\ ' 
Nj 


e 


/4tt 2 /3 
V N 2 


°)J 



(9) 



where the 9-function has the form 0(x, g) = ex P {~ ^J 2 + 27rigj}. Note that 

if the charge of the test particle is completely screened by the charge of the Higgs 
condensate (M/N 6 then the string tension is equal to zero, since Q(x, q + n) = 
Q(x,q) (n is an integer). The area law of the Wilson loop for fractionally charged 
particles in the 2D Abelian Higgs model was first found in Jj]], but then it was not 
realized that the nonzero value of the string tension is due to the analogue of the 
Aharonov-Bohm effect. 

For D = 3, 4, the expression (^) can be estimated in the saddle-point approxima- 
tion ||. Let us represent the condition S*j = in (|7]) by introducing integration over 
an additional field C: 



E 



*i(*c 2 )G^ *j(*c 2 )e^ 

S*j=0 



-\-oo 

E 5(5*j)---= J V*C £ ex V {i(5*j,*C)}- 



(10) 



oo *j(*c 2 )e^ 



^oo 

Inserting unity 1 = / V*F S (*F — *j) into equation (0), we get 

— oo 

+oo +oo 

< W M (C) >= J V*C J V*F ' 

-oo -oo *i(*C2)G^ 

exp \\*Ff + 2m— {*m c , *F) + *(5*F, *C)\-6 (*F - *j) . 



The use of the Poisson formula Ylm^z, S(n — x) = J22imez e%nx an d integration over the 
fields *F leads us to the dual representation of the quantum average (0): 



l r /V 2 M 

<Wm(C)>=zbkt J V *C £ ex P {-— 11^^ + 2^-^ + 2^-11 }. (12) 

-oo *i(*C2)e^ 



For 



iV 2 A (D 1 ) (0)»4/3, 



(13) 



where A^(R) is the D-dimensional massless lattice propagator, we can evaluate this 
expression in the semiclassical approximation: the result is 
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< W M (C) > N = canst. exp{-^(M, iV; /?) • (j c , ^jc)} , (14) 

where 



«(°)(M,JV;/3) = ±— -, g=min|— (15) 
v y 4/3 Rex 1 N 

Here g is the distance between the ratio M/N and the nearest integer number. Just as 
equation (§), the expressions flT4T ), (|15|) depend on the fractional part of M/N, which 
is a consequence of the Aharonov-Bohm effect. The interaction of the test charges is 
absent if q — (M/N is integer); this corresponds to the complete screening of the 
test charge M by the Higgs bosons of charge N. 

To find the potential induced by the Aharonov-Bohm effect we consider the product 
of two Polyakov lines: W M {C) = L+ (0) • L M (R). Then (j e , A _1 j c ) = 2T (A^_ 1} {R) + 

A7 D _ 1 -)(0)) and equation ( |T4T ) is reduced to 

<L+(0)L A/ (i?) >7V = const. exp{-2/ t (°)(M,iV;/?)TA^_ 1) ( J R)}. (16) 

For large R, we have A^(R) = f- lni?+. . ., and A^(R) = ^■R~ 1 + . . ., where c 3 and c 4 
are constants. Thus, the Aharonov-Bohm effect in three- and four-dimensional Abelian 
Higgs model gives rise in the leading order to the following long-range potentials: 

V^S(R) = 2 k (0) (M, N;/3) T A^ (R) = c 3 « (0) (M, N; (3) ■ In R + O^R' 1 ) , (17) 
V^iR) = 2 «< >(M, AT; P)TA^(R) = c 4 k (0) (M, AT; /?) • 1 + O (iT 2 ) . (18) 

A 3D vortex model in the continuum is discussed in Q, where it is shown that 
the potential has the form V£ D (R) ~ const. q 2 ip ln-j^, with ip being proportional to 
the vortex condensate, and R$ having the order of the vortex width. This result is in 
agreement with our estimation ( |T7f) of the potential. 

Note that from expression (|7|) the following properties of the potential V(m,n) can 
be obtained: 

V(M,N) = V(N-M,N) , V(N,N) = 0. (19) 

The long-range potentials V(R), equations flT7D, (0), and the string tension (0) satisfy 
these relations. 
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3 Numerical Calculations 



We calculated numerically the potential between the test particles of charge M in the 
three-dimensional Abelian Higgs model, the Higgs boson charge being N = 6. The 
action of the model is chosen in the Wilson form: S[A, (p] = /3\\dA\\ 2 — 7cos(dy? + 
NA). We perform the calculations by the standard Monte-Carlo method and we use 
200 statistically independent configurations for each point in the (3 — 7 plain. The 
simulations are performed on lattice of size 16 3 for the test charges M — 1, . . . , N. 

We fit the numerical data for the Wilson loop quantum average by the following 
formula: 



- In < W M (C) > N = K num (j c ,A 1 Jc ) + m 2 num \\j c \\ 2 + C num , (20) 

where n num , m num and C num are the fitting parameters. It turns out that the expec- 
tation values of the Wilson loops of the shapes L\ x L 2 , L%, L 2 = 1, . . . , 5 are well 
described by (f20|). 




I * • • • • 1 l l 1 l 1 l l 1 l 1 1 l 1 l 1 1 l 1 l 1 1 l 

0.8 1.2 1.6 2.0 2.4 2.8 

P 



Fig. 1: The ratio ( = ^^f vs. (3 for M = 1, 2, 3 at 7 = 1. 



In Fig. 1 we show the dependence of the coefficients ( = K num / on (3 for M = 
1,2,3 and 7 = 1, where is given by the semiclassical expression ([15|). We studied 
the region 1 < (3 < 3 for which the condition ( |13D of validity of the semiclassical 
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expression (|T5f) is not validP] (note that Aj^x(O) ~ 0.24.). The considered values of 
the parameters (3 and 7 correspond to the Coulomb phase in which the vortices are 
condensed. 

It is clearly seen from Fig. 1 that for (3 > 2.3 the values of ( are independent 
of M. Thus, n num is proportional to k^ ' with a coefficient which differs from unity 
because of quantum corrections. The same effect is observed in this model at finite 



temperature [10 



The numerical data also shows that = an d the potential for the charge 

M = N is equal to zero within numerical errors for all considered values of (3. These 
relations are in agreement with the Aharonov-Bohm nature of the potential V(m,n){R) 
(cf. equation (|19|)). Note that the usual Coulomb interaction of the test particles of 
charge M is proportional to M 2 . Therefore, the fact that the potential Vm) satisfies 
the relations means that the usual Coulomb interaction of the test particles is 
small. 



Conclusion and Acknowledgments 

We have shown by analytical calculations that the Aharonov-Bohm effect induces the 
long-range interaction (|^,|5]) at small values of (3: (3 < N 2 j, (3 < A (3 )(0)A^ 2 /4. Our 
numerical simulations for (3 > 2.3 show that this effect leads to long-range potential 
which is proportional to potential flT7[): V 3D (R; M, N, (3) = ((f3)V$ (R; M, N, (3), where 
( is independent of M. At the intermediate values of f3 the Aharonov-Bohm effect 
also induces a long-range potential but the induced potential is not proportional to 
the semiclassical expression eq. flTT|) . The potential depends non-analytically on the 
charge of the test particle. Due to the long-range nature of the induced potential, the 
Aharonov-Bohm effect may be important for the dynamics of the colour confinement 
in nonabelian gauge theories |TT| . 

M.N.Ch. and M.I. P. acknowledge the kind hospitality of the Theoretical Depart- 
ment of the Kanazawa University. The authors are grateful to D.A. Ozerov for useful 
discussions. This work is supported by the JSPS Program on Japan - FSU scientists 
collaboration, and also by the Grants: INTAS-94-0840, INTAS-94-2851, INTAS-RFBR- 
95-0681, and Grant No. 96-02-17230a of the Russian Foundation for Fundamental 
Sciences. 



2 Condition (|T^) for N — 6 reads as follows: (3 <C 2.2. Unfortunately, even at j3 < 1 our data for £ 
contains large error bars, therefore we do not show this region of f3 on the figure. 
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Appendix A 



Let us briefly summarize the main notions from the theory of differential forms on the 
lattice [fj. The advantage of the calculus of differential forms consists in the general 
character of the expressions obtained. Most of the transformations depend neither on 
the space-time dimension nor on the rank of the fields. With minor modifications, the 
transformations are valid for lattices of any form (triangular, hypercubic, random, etc). 
A differential form of rank k on the lattice is a function defined on k- dimensional 
cells Cfc of the lattice, e.g., the scalar (gauge) field is a 0-form (1-form). The exterior 
differential operator d is defined as follows: 

(<ty)(C fc+ i) = E Wu). (A.l) 

C k &9C k+l 

Here dCk is the oriented boundary of the fc-cell C^. Thus, the operator d increases the 
rank of the form by unity; dip is the link variable constructed, as usual, in terms of 
the site angles (p, and dA is the plaquette variable constructed from the link variables 
A. The scalar product is defined in the standard way: if (p and ip are fc-forms, then 
(<f,ip) = I2c k ( p(C'k)'4 ) (C'k)i where J2c k * s the sum over all cells C^. To any /c-form 
on the .D-dimensional lattice there corresponds a (D — fc)-form *$(*Cfc) on the dual 
lattice, *Ck being the (D — fc) -dimensional cell on the dual lattice. The co-differential 
5 = *d* satisfies the partial integration rule ((p,5ip) = (d(p,ip). Note that 5$(Cfc) is 
a (k — l)-form and 5$(Co) = 0. The norm is defined by ||a|| 2 = (a, a). Therefore, 
\\dtp + 2ttI — NA\\ 2 in (^j implies summation over all links. The sum 52i(ci)ez ^ s taken 
over all configurations of the integers I attached to the links C\. Due to the well-known 
property d 2 = 5 2 = 0, the action (Q) is invariant under the gauge transformations 
A' = A + da, <p' = ip + Na . The lattice Laplacian is defined by A = d5 + 5d. 
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